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Large  cylindrical  test  beds  2.13  m  (7  ft)  in  diameter  and  cast  from  SRl- 
RMG-2C4  grout  were  used  by  Merritt  CASES,  Inc.,  and  Physics  International  Com¬ 
pany  to  investigate  the  simulation  of  nuclear  test  effects  on  deep-based 
structures  in  the  "High-Explosive  Model  Structure  Simulation  (HEMSS)  Program." 
2C4  grout  was  used  in  the  test  beds  and  reached  temperatures  in  excess  of  108°C 
(226"F)  in  HEMSS  No.  1  and  2  cylinders.  Although  the  grout  properties  were 
well  documented,  the  effects  of  heat  of  hydration  on  the  grout  (Continued) 
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,0.  ABSTRACT  (Continued) 

when  cast  in  large  volumes  were  unknown. 

Thermal  studies  of  tlie  llEMSS  cylinders  were  conducted  at  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  (WES)  based  upon  computer  simulations. 
These  studies  were  conducted  to  determine  the  tliermal  response  of  2C4  grout 
placed  in  mass._.  The  computer  simulations  were  conducted  in  four  groups.  Group 
I  runs  were  conducted  to  duplicate  the  thermal  histories  in  the  HEMSS  No.  1 
cylinder,  including  the  establishment  of  tliermal  properties.  Group  II  computer 
runs  were  made  to  evaluate  several  temperature  control  measures  which  could  be 
employed  on  future  HEMSS  test  beds  to  reduce  internal  temperature  and  potential 
excessive  thermal  stress  or  strain.  Measures  recommended  included  reducing 
placement  temperature  and  ventilating  the  0.3048-m  (1-ft)  hole  at  the  axis  of 
the  cylinders.  Group  111  runs  predicted  temperatures  in  the  HEMSS  No.  3  cylin¬ 
der  while  evaluating  the  effects  of  the  diurnal  ambient  temperature  cycle  and 
the  several  attainable  placement  temperatures.  Group  IV  runs  include,  postcast¬ 
ing  temperature  and  thermal  stress-strain  evaluation  of  the  2.13-m  (7-ft)  diam¬ 
eter  HEMSS  No.  3  cylinder  as  well  as  the  same  analysis  on  potential  3.66-  and 
6.10-m  (12-  and  20-ft)  diameter  cylinders. 


The  study  showed  that  the  finite-element  method  temperature  simulation 
program  in  use  at  WES  could  predict  temperature  distribution  in  the  HEMSS  cylin¬ 
ders  as  a  function  of  time  after  casting.  Temperature  in  the  2.13-m  (7-ft) 
diameter  cylinders  could  only  be  effectively  controlled  by  lowering  the  place¬ 
ment  temperature  of  the  grout.  Simulations  of  3.66-  and  6.10-m  (12-  and  20-ft) 
diameter  cylinders  indicated  increasing  peak  temperature  and  thermal  strains  as 
cylinder  diameter  increased.  With  conditions  similar  to  those  of  the  HEMSS  No. 

3  cylinder,  a  6.10-m  (20-ft)  diameter  cylinder  is  predicted  to  reach  a  peak 
temperature  of  116°C  (240'’F).  It  is  recommended  that  additional  work  be  done 
to  document  the  properties  of  2C4  grout  when  cured  at  temperatures  at  or  near 
100°C  (212‘’F). 
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This  study  was  conducted  in  the  Structures  Laboratory  (SL)  of  the 
U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  under  the  sponsor¬ 
ship  of  the  Defense  Nuclear  Agency  as  part  of  the  FY  1978  Nuclear  Weapons 
Effects  Program,  Subtask  J34HAXSX311,  "Underground  Structures." 

The  study  was  conducted  under  the  general  supervision  of  Mr. 

Bryant  Mather,  Chief,  SL,  and  Mrs.  Katharine  Mather,  former  Chief,  Engi¬ 
neering  Sciences  Division,  and  under  the  direct  supervision  of  Mr.  Billy 
R.  Sullivan,  former  Chief,  Engineering  Physics  Branch.  Mr.  Anthony  A. 
Bombich  was  project  leader  and  prepared  this  report. 

Funds  for  publication  of  the  report  were  provided  from  those  made 
available  for  operation  of  the  Concrete  Technology  Information  Analysis 
Center  (CTIAC) .  This  is  CTIAC  Report  No.  56. 

Commander  and  Director  of  WES  when  this  report  was  published  was 
COL  Tilford  C.  Creel,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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Conversion  Factors,  i{on-SI  to  SI  (Metric) 

Units  of  Measurement 

Non-Sl  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply _  _ ^ _  _ To  Obtain 


British  thermal  unit 
(International  table)  •  feet 
per  hour  •  square  foot  • 
degree  Fahrenheit 

1.730735 

watts  per  metre  • 
Kelvin 

British  thermal  units 
(International  table)  per 
pound  (mass)  *  degree 
Fahrenheit 

4186.8 

joules  per 
kilogram  •  Kelvin 

calories  (International  table) 

4186.8 

joules  per  kilogram 

per  gram 

cubic  feet  per  minute 

0.0004719474 

cubic  metres  per 
second 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or 
Kelvins* 

feet 

0.3048 

metres 

feet  per  second 

0.3046 

metres  per  second 

inches 

0.0254 

metres 

inches  per  inch  per  degree 
Fahrenheit 

5/9 

metres  per  metre 
per  Kelvin 

miles  (U.  S.  statute)  per  hour 

1.609347 

kilometres  per  hour 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (force)  per  square  inch 

6894.757 

pascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic 
metre 

square  feet  per  hour 

0.0000258064 

square  metres  per 

second 


*  To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F) 
readings,  use  the  following  formula;  C  =  (5/9)(F  -  32).  To  obtain 
Kelvin  (K)  readings,  use:  K  =  (5/9)(F  -  32)  +  273.15. 
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THERMAL  STUDIES  OF  HEMSS 


CYLINDRICAL  TEST  BEDS 
CONTAINING  SRI-RMG-2C4  GROUT 


Introduction 


1.  This  study  was  conducted  at  the  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  to  evaluate  the  heat  generated  by  SRI-MRG-2C4 
grout  cast  in  large  cylindrical  test  beds.  Merritt  CASES,  Inc.,  and 
Physics  International  Company  are  conducting  the  "High-Explosive  Model 
Structure  Simulation  (HEMSS)  Program"  for  the  Defense  Nuclear  Agency 
(DNA)  to  investigate  the  simulation  of  nuclear  test  effects  on  deep- 
based  structures  using  high  explosives  (HE)  in  rock  simulant  test  beds. 
The  2C4  grout  is  a  tuff  simulant.  This  grout,  originally  developed  at 
WES  and  later  modified  by  Stanford  Research  Institute  (SRI) ,  was  chosen 
from  among  several  other  grout  mixtures  after  comparing  its  mechanical 
properties  with  those  of  Nevada  Test  Site  tuff.  Its  static  response 
is  well  documented.  However,  the  effects  of  heat  of  hydration  when  the 
grout  is  cast  in  large  volumes  were  unknown.  The  composition  in  percent 


by  weight  of  SRI-RMG-2C4  is: 

Type  G  portland  cement  (or  Type  l-II)  32.691 
Sand  (20-40)  Monterey  21.896 
Barite  (barium  sulfate)  20.848 
Bentonite  (gel)  2.837 
D65  (concrete  friction  reducing  compound)  0.078 
Water  21.650 


2.  When  this  study  began,  two  cylinders  designated  HEMSS  No.  1 
and  No.  2,  2.13  m  (7  ft)*  in  diameter  and  4.27  m  (14  ft)  in  height,  had 
been  cast.  A  0.3048-m  (1-ft)  diameter  vertical  hole  was  located  concen¬ 
tric  with  the  axis  of  each  cylinder.  It  was  understood  that  the  holes 
were  not  vented.  The  assumption  was  made  at  WES  that  near-adiabatic 

*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  preser.ted  on  page  3. 
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conditions  existed  at  the  hole  at  the  inner  surface  of  the  cylinders. 
Thermistors  had  been  used  to  measure  temperatures  in  the  cylinders  for 
approximately  the  first  6  days  after  placement.  The  thermistors  were 
located  on  a  radius  at  equal  increments  of  0.229  m  (9  in.)  between  the 
inner  and  outer  surfaces  of  the  cylinders.  Peak  temperatures  recorded 
in  HEMSS  No.  1  and  No.  2  at  a  distance  of  0.229  m  (9  in.)  from  the 
inner  surface  were  108  and  128  C  (226  and  262  F) . 

3.  The  objectives  of  this  study,  narrowed  from  those  originally 
proposed  because  of  a  reduced  level  of  funding,  included  efforts  to 
duplicate  analytically  the  thermal  histories  of  the  HEMSS  cylinders. 

In  subsequent  simulations  it  was  intended  to  determine  the  effects  of 
several  temperature  control  measures  that  could  be  used  in  subsequent 
castings.  Finally  it  was  intended  to  conduct  simulations  of  temperature 
history  and  resulting  thermal  strain  in  the  2,13-m  (7-ft)  diameter  and 
in  proposed  3.66-  and  6.10-m  (12-  and  20-ft)  diameter  cylinders, 

4.  Two  computer  programs  were  used  in  the  study.  These  are 

1  2 

two-dimensional  finite-element  programs  ’  developed  at  the  University 
of  California,  Berkeley,  by  Wilson,  Sandhu,  and  associates  for  the 
U,  S.  Army  Engineer  District,  Walla  Walla.  These  programs  were  developed 
to  compute  thermal  histories  and  resulting  stress-strain  histories  within 
mass  concrete  structures.  The  programs  permit  simulation  of  incremental 
construction  (placement  of  concrete  by  lifts),  allow  varying  boundary 
conditions,  and  require  time-dependent  heat  generation  and  mechanical 
properties . 

5.  Computer  simulations  for  HEMSS  cylinders  were  conducted  in 
four  groups  of  runs. 

Group  I.  Duplicate  the  results  of  HEliSS  No.  1  including 
establishing  estimated  thermal  properties. 

b.  Group  II.  Evaluate  the  effects  of  several  temperature 
control  measures. 

£.  Group  III.  Temperature  predictions  for  HEMSS  No.  3. 

Group  IV.  Temperature  and  thermal  stress-strain  analysis 
of  2.13-,  3.66-,  and  6.10-m  (7-,  12-,  and  20-ft)  diameter 
cylinders  based  upon  conditions  obtained  from  HEMSS  No.  3 
measurements . 
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Temperature  Calculations 

6.  It  was  only  feasible  to  do  a  limited  amount  of  testing  to 
establish  the  thermal  properties  of  2C4  grout  for  computer  input.  It 
was  therefore  necessary  to  estimate  the  properties  based  upon  information 
available  such  as  the  mixture  proportions  and  materials  and  measured 
temperatures  in  the  HEMSS  No.  1  cylinder.  Very  little  documentation  of 
thermal  properties  exists  for  grouts.  Estimates  of  2C4  thermal  properties 
were  modifications  of  known  tliermal  behavior  of  concrete.  It  is  known 
that  the  specific  heat  of  concrete  increases  as  temperature  increases. 

It  was  suspected,  and  while  the  report  was  in  progress,  learned  that  the 
specific  heat  of  the  cement  and  water  undergoes  an  irreversible  change 
during  hydration.  The  work  documenting  this  phenomenon  has  not  yet  been 
received;  however,  it  is  known  that  the  specific  heat  is  lowered  during 
hydration.  In  conventional  concrete  the  total  effect  of  the  change  is 
very  small  because  tt)e  percentage  of  water  by  mass  of  the  mixture  is  very 
low.  This  change  in  effect  reduces  the  specific  lieat  of  water  from 
4186. 8  J/kg-K  (1.0  Btu/lb-“F)  to  837.4-1256.0  J/kg-K  (0.2-0. 3  Btu/lb-°F) 
or  near  to  the  specific  heat  of  tlie  cement.  In  grouts,  and  in  particular 
in  2C4  grout,  the  water  content  by  mass  is  greater  than  20  percent. 

Since  the  cement  content  is  very  high,  the  effect  on  reducing  the  specific 
heat  of  the  mixture  due  to  the  hydration  phase  change  is  greater  than  in 
concrete. 

7.  The  specific  heat  of  concrete  or  grouts  is  also  highly  tempera- 

3 

ture  dependent.  For  example,  specific  heat  of  concrete  increases  approxi¬ 
mately  20  percent  as  the  temperature  Increases  from  21  C  (70  F)  to  66  C 
(150  F)  with  approximately  20  percent  additional  increase  to  93  C  (200  F) . 
Although  the  increase  is  loss  in  mixtures  containing  more  water,  as  is 
the  case  with  grout,  one  may  assume  a  30-40  percent  increase  in  specific 
heat  for  2C4  grout  between  placement  at  21  C  (70  F)  and  peak  temperatures 
over  93  C  (200  F) .  The  implications  of  temperature  dependence  and  hydra¬ 
tion  phase  change  of  specific  lieat  in  this  study  are  large.  During  the 
initial  period  when  intense  hydration  is  occurring,  the  specific  heat  is 
increasing  due  to  the  temperatu’‘e  rise  and  decreasing  with  Increasing 
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hydration  or  age.  As  a  result  of  these  two  inverse  functions  of  tempera¬ 
ture  and  age  during  temperature  rise,  the  heat  required  to  raise  the 
temperature  will  be  greater  than  the  heat  lost  as  the  temperature  drops 
an  equivalent  amount, 

8.  Prior  studies  at  WES  which  included  computer  calculations  of 
temperature  histories  in  structures  containing  portland  cement  concrete 
mixtures  have  not  considered  temperature-dependent  thermal  properties 
because  the  lean  concrete  mixtures  used  and  substantial  temperature 
controls  employed  do  control  temperature  rise  and  temperature  dependence. 
For  this  reason  computer  programs  currently  used  were  prepared  to  con¬ 
sider  only  age-dependent  properties.  The  temperature-  and  age-dependent 
properties  characteristic  of  the  materials  and  conditions  in  the  HEMSS 
cylinders  could  not  be  incorporated  without  additional  computer  program 
development.  Therefore,  thermal  properties  were  established  for  2C4 
grout  which  would  yield  acceptable  comparative  thermal  simulation  with 
the  HEMSS  cylinders,  A  brief  discussion  of  the  development  of  input 
thermal  properties  to  the  computer  program  follows.  A  value  of  1457  J/ 
kg’K  (0.348  Btu/lb*°F)  for  specific  heat  was  used.  This  is  the  initial 
theoretical  value  calculated  from  the  mixture  proportions  and  the  specific 
heat  of  each  component  material.  A  value  for  thermal  conductivity  was 
obtained  as  follows.  Conductivity  is  the  product  of  specific  heat, 
thermal  diffuslvity,  and  density.  Thermal  diffuslvity  was  determined 

in  accordance  with  CRD-C  36.^  Three  76.2-  by  152.4-mm  (3-  by  6-in.) 

cylinders  of  2C4  grout  were  cast  with  copper  constantan  thermocouples 

embedded  at  the  centroid  of  each  cylinder.  Diffuslvity  conducted  on 

-7  2  2 

the  cylinders  at  3-days  age  yielded  4.129  x  10  m  /sec  (0.016  ft  /hr). 

3  3 

A  density  of  2160  kg/m  (134.78  Ib/ft  )  was  chosen  from  data  provided 
by  CASES.  Calculated  thermal  conductivity  was  1.305  W/m-K  (0.754  Btu* 
ft/hr-ft^-“F). 

9.  The  additional  thermal  property  required  is  adiabatic  tempera¬ 
ture  rise  which  provides  the  heat  of  hydration  history.  Experimental 
testing  for  adiabatic  temperature  rise  was  not  Included  originally 
because  it  was  known  that  the  high  temperatures  of  2C4  grout  exceeded 
the  operating  limit  of  existing  facilities.  During  the  course  of  the 
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study  other  existing  equipment  was  being  modified  to  allow  at  least 
abbreviated  adiabatic  tests  of  mixtures  of  high  temperature  rise.  Satis¬ 
factory  operations  were  not  achieved  in  time  for  this  study.  An  adiabatic 
temperature  rise  history  was  then  obtained  from  the  measured  temperature 
in  HEMSS  No.  1.  It  was  assumed  that  temperatures  measured  during  the 
first  day  after  placement  at  the  location  nearest  to  the  axis  of  the 
cylinder  were  adiabatic.  To  establish  the  adiabatic  history  after  the 
first  day  it  was  necessary  to  establish  the  heat  of  hydration  of  the 
first  day  relative  to  the  total  heat  generation  possible.  Using  the 
measured  temperature  rise  of  66  C  (150  F)  after  the  first  day,  a  specific 
heat  of  1457  J/kg-K  (0.345  Btu/lb"“F)  for  the  grout,  and  the  specified 
cement  content,  the  heat  of  hydration  calculated  was  0.364  Mj/kg  (87  cal/g). 
This  is  greater  than  that  typical  for  28-day  Isothermal  heat  of  hydration 
test  results  of  Type  "G”  cement;  it  is  also  approximately  twice  that  nor¬ 
mally  measured  at  1-day  age.  Typical  28-day  heat  of  hydration  values  are 
perhaps  80-90  percent  of  the  potential  heat  generation  possible  in  stand¬ 
ard  tests.  It  was  concluded  that  because  of  the  initially  high  tempera¬ 
tures  accelerating  the  rate  of  hydration  that  only  an  additional  10  per¬ 
cent  of  the  total  heat  would  be  generated  by  28  days  and  that  by  28  days 
hydration  would  have  ceased.  The  remaining  adiabatic  temperature  rise 
history  after  first  day  was  then  fabricated  assuming  an  11  C  (20  F)  rise 
between  1  and  28  days  age. 

Group  I  temperature  runs 

10.  Runs  1-6  were  made  to  establish  that  input  data  chosen  were 
sufficiently  correct  to  duplicate  the  measured  temperatures  in  HEMSS  No.  1. 
A  finite-element  model  was  prepared.  Tliis  model  (Figure  1)  was  a  segment 
of  a  horizontal  plane  tlirough  tlie  cylinder.  Heat  was  assumed  to  flow 
radially  or  only  througli  tlie  inner  or  outer  cylindrical  surfaces.  The 
vertical  hole  in  tiie  center  of  the  cylinder  was  assumed  to  be  sealed  from 
outside  air;  thus,  in  runs  1-6  heat  was  assumed  lost  only  through  the 
outer  surface.  In  runs  1-6  input  properties  and  ambient  air  /elocities 
used  for  the  convection  heat-transfer  coefficient  at  the  outside  surface 
were  varied  slightly  to  obtain  a  set  of  input  data  that  produced  best 
correlation  with  measured  t emper  itures .  Placement  temperature  was 


8 


assumed  to  be  24  C  (75  F)  and  ambient  air  was  assumed  still  and  23  C 
(73  F) .  Ambient  air  velocity  was  ultimately  increased  to  0.0762  m/sec 
(2.5  ft /sec)  in  run  6.  Table  1  lists  the  pertinent  input  data  for  runs 
1-6  as  well  as  all  other  runs  in  this  study.  Figure  2  shows  calculated 
temperatures  in  run  6  versus  measured  temperatures  in  HEMSS  No.  1.  The 
slower  rate  of  temperature  drop  in  run  6  when  compared  to  HEMSS  No.  1 
after  2  days  probably  results  from  some  or  all  of  the  following: 
incorrect  assumed  air  velocity,  use  of  constant  rather  than  temperature 
dependent  thermal  properties  in  the  calculations,  and  heat  losses  to  the 
ends  of  the  cylinder  where  no  such  loss  was  assumed  in  the  calculations. 
Input  data  used  in  run  6  were,  however,  selected  as  the  best  set  for  sub¬ 
sequent  runs.  Adiabatic  temperature  rise  used  in  run  6  as  well  as  all 
subsequent  runs  is  shown  in  Figure  3. 

Group  II  temperature  runs 

11.  Runs  7-9  were  made  to  evaluate  several  potential  temperature 
control  measures.  One  method  was  to  reduce  placement  temperature.  Run  7 
was  similar  to  run  6  except  that  placement  temperature  was  reduced  from 
24  C  (75  F)  to  10  C  (50  F) ,  Adiabatic  conditions  were  assumed  at  the 
inner  surface.  Run  7  showed  that  a  14  C  (25  F)  reduction  In  placement 
temperature  reduced  peak  temperatures  by  approximately  the  same  amount, 
or  108  C  (226  F)  rec’uced  to  94  C  (202  F) .  Past  experience  has  shown  that 
an  equivalent  reduction  between  peak  and  placement  temperatures  will 
occur  only  when  peak  temperature  is  reached  very  soon  after  placement  or 
when  substantial  time  passes  between  placing  sequential  lifts  in  a 
structure.  Run  8  was  conducted  to  examine  the  effects  of  forced  ventila¬ 
tion  on  the  surface  of  the  0.3048-m  (1-ft)  diameter  hole.  It  is  obvious 
because  of  geometry  that  given  equal  air  temperature  and  air  velocity 
that  substantially  more  heat  will  flow  from  the  outside  than  from  the 
inside  surface  of  a  hollow  cylinder.  The  smaller  the  ratio  of  radii  of 
the  inner  and  outer  surfaces,  the  greater  is  this  difference.  Since  the 
cylinder  rests  on  end,  ventilation  of  tiie  hole  can  only  occur  from  the 
top  surface.  In  run  8  a  ventilation  pipe  was  assumed  positioned  in  the 

hole;  tlie  bottom  of  the  pipe  was  assumed  a  foot  or  so  above  the  bottom 
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of  the  hole.  A  blower  having  a  capacity  of  0.472  m  /sec  (1000  ft  /min) 


was 


assumed  to  blow  down  the  pipe  venting  upward  between  the  pipe  and 
surface  of  the  hole.  A  pipe  diameter  of  0.203  m  (8  in.)  was  chosen 
because  the  cross-section  for  inlet  and  outlet  were  equal  in  a  0.3048-m 
(12-in.)  diameter  hole.  Based  upon  these  dimensions  and  the  air  flow 
specified,  a  velocity  of  13.4  m/sec  (44  fps)  was  assumed  to  ventilate 
the  inner  surface  of  the  cylinder.  All  other  parameters  were  the  same 
as  run  6.  Peak  temperature  in  run  8  was  102  C  (215  F)  occurring 
equidistant  between  the  inner  and  outer  surface  of  the  cylinder  at  1,0 
day.  This  is  0.25  days  earlier  and  1.7  C  (3  F)  less  than  occurred  in 
run  6  at  the  same  location.  Comparison  of  temperatures  near  the  inner 
surface  of  the  cylinders  shows  108  C  (226  F)  in  run  6  at  1.25-days  age 
versus  37  C  (99  F)  occurring  at  0.5-days  age  in  run  8  during  forced 
ventilation.  Run  8  demonstrates  a  substantial  reduction  in  peak  tempera¬ 
ture  at  the  inner  surface  even  though  peak  temperature  within  the  cylinder 
is  not  lowered  very  much.  The  position  at  which  peak  temperature  occurs 
may  be  important.  Although  a  thermal  stress  simulation  was  not  made  in 
this  study  for  the  unvented  hole,  it  is  probable  that  when  highest 
temperatures  occur  at  the  center  of  the  cylinder,  thermal  stress  will  be 
higher  than  in  one  cooled  at  the  center, 

12.  Run  9  was  made  using  the  10  C  (50  F)  placement  temperature 
with  forced  ventilation  at  the  inner  surface.  Peak  temperature  was  89  C 
(192  F)  occurring  at  the  same  location  as  in  run  8.  The  reduction  in 
peak  temperature  was  equivalent  to  the  reduction  in  placement  temperature. 
The  results  of  runs  7-9  (Figures  4-6)  show  that  a  reduction  in  placement 
temperature  is  effective  in  reducing  peak  temperatures.  Although  forced 
ventilation  at  the  inner  surface  does  not  substantially  reduce  tempera¬ 
tures  at  a  point  equidistant  between  tlie  inner  and  outer  surfaces,  tem¬ 
peratures  at  the  inner  surface  are  reduced  more  than  50  percent.  It  was 
recommended  to  Merritt  CASES,  Inc.,  that  both  control  measures  be  employed 
on  future  HEMSS  cylinders. 

Group  III  temperature  runs 

13.  These  runs  were  requested  by  Merritt  CASES,  Inc.,  in  support 
of  casting  the  HEMSS  No.  3  cylinder.  The  objectives  were  to  evaluate  the 
effects  of  the  diurnal  ambient  temperature  cycle,  different  times  of 
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placement  during  the  day,  and  the  several  attainable  placement  tempera¬ 
tures.  HEMSS  No.  3  was  placed  8  August  1978.  Estimates  were  made  by 
Merritt  CASES,  Inc.,  for  realistic  placement  temperatures  attainable  in 
August  using  chilled  mixing  water.  These  were  13  C  (55  F)  and  18  C 
(65  F)  for  placements  beginning  at  6  a.m.  and  6  p.m. ,  respectively.  The 
higher  afternoon  temperatures  are  due  to  solar  heating  of  the  materials 
and  equipment  during  the  day.  It  was  assumed  that  placements  of  the 
grout  at  midheight  of  the  cylinder  would  occur  3  hr  after  placement 
began.  This  is  the  assumed  location  of  the  temperature  simulations. 

Figure  7  shows  the  estimated  diurnal  ambient  temperatures  provided  by 
Merritt  CASES,  Inc.,  that  were  used  in  runs  10-12.  Runs  10  and  11  were 
simulations  of  morning  and  evening  placement.  Zero  time  on  the  ambient 
temperature  curves  were  9  a.m.  and  9  p.m. ,  respectively,  and  placement 
temperatures  were  13  C  (55  F)  and  18  C  (65  F) ,  respectively.  The  inner 

Surface  was  ventilated  with  ambient  air  at  9.1  m/sec  (30  ft/sec),  an 
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estimate  of  actual  air  velocity  from  a  0.425  m  /sec  (900  ft  /min)  blower 
assuming  some  losses  due  to  the  180-degree  turn  at  the  bottom  of  the  hole. 
Air  velocity  assumed  at  the  outer  surface  of  the  cylinder  was  4.48  m/sec 
(10  raph) ,  which  is  the  daily  average  velocity.  Solar  radiation  effects 
were  not  considered. 

14.  Peak  temperatures  were  92  C  (198  F)  and  97  C  (206  F)  for 
runs  10  and  11,  respectively,  midway  between  the  inner  and  outer  surfaces, 
or  1.1  C  (2  F)  less  than  the  difference  in  placement  temperatures. 
Observing  that  these  differences  are  nearly  equivalent,  it  can  be  con¬ 
cluded  that  placement  during  the  morning  subjecting  tlie  mass  initially 
to  the  higher  temperatures  during  the  afternoon  does  not  affect  peak 
temperatures  when  compared  to  evening  placement  followed  by  cooler  night¬ 
time  temperatures.  This  applies  to  the  case  when  direct  solar  radiation 
is  Ignored.  Therefore,  the  criterion  that  should  govern  selection  of 
the  best  placement  time  is  the  one  that  most  easily  allows  lowest  place¬ 
ment  temperatures.  The  effects  of  direct  solar  radiation  on  the  outer 
surface  steel  forms  could  change  this  conclusion  if  the  forms  are  suffi¬ 
ciently  heated  to  reduce  early  cooling.  However,  this  is  not  likely  for 
2C4  grout  cast  in  at  least  the  dimensions  of  HEMSS  No.  1-No.  3.  The  rate 
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of  initial  hydration  overpowers  the  boundary  condition. 

15.  Run  12  was  the  same  as  run  10;  that  is,  morning  placement  at 
12.8  C  (55  F)  except  that  adiabatic  conditions  were  maintained  at  the 
inner  surface.  Peak  temperature  at  the  midpoint  between  inner  and  outer 
surfaces  was  94  C  (202  F)  or  2.2  C  (4  F)  higher  than  in  run  10  at  the 
same  location.  However,  inner  surface  temperatures  were  97  C  (207  F) 
compared  with  47  C  (117  F)  in  run  10.  Although  peak  temperatures  differ 
by  only  5  C  (9  F)  between  runs  10-12,  positive  ventilation  should  be 
used  because  less  volume  in  the  cylinders  is  subjected  to  the  excessively 
high  temperatures  and  in  all  probability  potential  thermal  stress  is  less. 
Time  histories  of  temperature  change  for  these  runs  are  shown  in  Figures 
8-10  at  the  six  thermistor  locations  specified  for  HEMSS  No.  3.  These 
data  were  discussed  by  telephone  with  Merritt  CASES,  Inc.,  including 
recommendations  for  casting  the  HEM.SS  No.  3  cylinder.  Pertinent  data 
from  runs  7-12  were  forwarded  to  Merritt  CASES,  Inc.,  on  15  September 
1978. 

Group  IV  temperature  runs 

16.  Following  receipt  of  temperatures  measured  in  HEMSS  No.  3, 
several  simulations  (runs  13-19)  were  made  to  duplicate  the  results,  to 
evaluate  the  boundary  conditions  to  include  solar  radiation,  and  to 
provide  the  basis  for  comparative  simulations  including  thermal  stress 
analysis  for  2.13-,  3.66-,  and  6.10-m  (7-,  12-,  and  20-ft)  diameter 
cylinders.  The  table  below  contains  tlie  radial  locations  of  thermistors 
provided  by  Merritt  CASES,  Inc.  These  thermistor  designations  will  be 
used  in  the  discussion  that  follows. 

Thermistor  _ Radius 

T1  0.15  m  (6.0  in.) (located  on  inner  form  in  air  flow) 

T2  0.23  m  (9.0  in.) 

T3  0.46  m  (18.0  in.) 

T4  0.61  m  (24.0  in.) 

T5  0.76  m  (30.0  in.) 

T6  1.07  m  (42.0  in.) (located  on  outer  steel  form) 

17.  The  measured  temperature  histories  in  HEMSS  No.  3  are  shown 
in  Figure  11.  Several  observations  were  possible.  First,  acknowledging 
that  the  surface  temperature  w<nild  not  be  correct,  run  13  was  made  to 
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see  if  by  correcting  run  10  for  the  actual  placement  temperature  of 
25.5  C  (78  F)  whether  peak  temperatures  would  be  duplicated.  Peak  tem¬ 
perature  in  run  13  was  102  C  (216  F)  comparing  with  101.7  C  (215  F)  in 
HEMSS  No.  3.  Secondly,  thermistor  T1  located  on  the  inner  surface  in 
the  air  flow  recorded  temperatures  approximately  8-11  C  (15-20  F)  below 
ambient  temperature  at  all  times.  The  question  that  arose  was  did  the 
thermistor  measure  air  temperature  or  is  the  temperature  measured  the 
effect  of  phenomena  that  lowered  the  thermistor  temperature  below  that 
of  the  air  such  as  might  be  expected  with  evaporative  cooling.  This 
question  gains  Importance  when  crying  to  establish  the  actual  boundary 
conditions  for  computer  simulation.  Runs  13-19  established  that  the  inner 
surface  air  temperature  was,  as  measured,  lower  than  ambient  air  tempera¬ 
ture.  Considerable  speculation  arose  as  to  the  cause;  however,  since  no 
cause  can  be  proven,  the  subject  was  dropped. 

18.  The  last  point  to  be  discussed  pertaining  to  the  HEMSS  No.  3 
temperatures  is  the  outer  surface  temperatures.  Gage  T6  was  located  on 
the  southern  face.  Therefore,  it  recorded  surface  temperatures  that 
Included  the  grout  temperature  plus  temperature  rise  due  to  solar  heating. 
After  examining  Figure  11  it  is  reasonable  to  conclude  that  the  surface 
temperature  history  exclusive  of  solar  effects  lies  along  a  line  of  best 
fit  through  the  minimum  temperatures  reached  at  night.  The  effects  due 
to  diurnal  air  temperature  is  a  temperature  swing  of  perhaps  1.6-5. 5  C 
(3-10  F)  above  this  line  at  approximately  2  p.m.  If  we  assume  an  air 
temperature  range  between  29  C  (85  F)  and  40  C  (105  F)  for  the  first 
several  days  after  placement,  followed  by  a  reported  cooling  resulting 
in  temperatures  between  18  C  (65  F)  and  29  C  (85  F)  for  several  additional 
days,  then  the  peak  surface  temperatures  are  slightly  below  the  peak  sol- 
air  temperatures.  Scl-air  temperatures  are  used  by  the  American  Society 
of  Heating,  Refrigeration  and  Air-Conditioning  Engineers  (ASHRAE).^  Sol- 
air  temperature  is  a  combination  of  ambient  temperature  and  solar  effects 
to  represent  an  equivalent  ambient  temperature  that  i.  used  to  calculate 
cooling  requirements  in  buildings.  The  use  of  the  sol-air  concept  allows 
the  convection  heat-transfer  method  in  the  temperature  program  to  account 
for  solar  radiation.  This  was  the  first  opportunity  at  WES  to  use  the 
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sol-air  concept  having  measured  effects  in  hand  to  verify  results, 

19.  Merritt  CASES,  Inc.,  requested  that  simulations  on  the  2,13-, 
3.66-,  and  6.10-m  (7-,  12-,  and  20-ft)  cylinder  be  made  using  the  actual 
conditions  deduced  for  HEMSS  No,  3,  Also  specified  to  be  used  was  a 
0.61-m  (2-ft)  diameter  axial  hole  in  the  3.66-m  and  6,10-m  (12-  and  20-ft) 
diameter  cylinders  exposed  to  the  same  temperature  environment  and  forced 
air  velocity  as  the  2.13-m  (7-ft)  cylinder  that  had  a  0.3048-m  (1-ft) 
diameter  hole.  All  earlier  temperature  runs  were  calculated  through 

7  days,  the  maximum  duration  of  measured  temperatures  after  placement. 
Since  the  2.13-m  (7-ft)  cylinder  did  not  completely  achieve  temperature 
equilibrium  during  this  time,  it  was  decided  that  the  2.13-m  and  3.66-m 
(7-  and  12-ft)  cylinder  simulations  should  be  run  for  14  days  and  the 
6.10-m  (20-ft)  cylinder  for  28  days.  The  boundary  ambient  conditions 
used  in  these  runs  were  the  modified  results  of  gages  T1  and  T6  checked 
out  in  trial  runs  13-19  to  7  days.  After  7  days,  constant  inner  and 
outer  ambient  temperatures  of  18  C  (65  F)  and  29  C  (85  F) ,  respectively, 
to  the  duration  of  each  run  were  used.  Forced  cooling  of  the  inner 
surface  was  assumed  to  continue  after  7  days. 

20.  A  finite-element  model  was  prepared  for  each  cylinder.  Each 
model  comprised  an  arc  of  90  deg  to  allow  displacement  restraint  to  be 
applied  in  stress  analysis  runs.  Each  model  consisted  of  the  same  number 
of  nodes  and  elements.  Nodes  were  located  radially  to  provide  the  best 
description  of  temperature  gradient  and  thermal  stress  and  strain. 

Figure  12  shows  the  model  for  the  2.13-m  (7-ft)  cylinder.  Figure  13 
shows  the  radial  node  and  element  centroid  layouts  used  for  the  three 
cylinder  diameters. 

21.  Run  20  was  the  final  simulation  of  HEMSS  No.  3.  Figures  14 
and  15  show  calculated  temperatures  on  two  time  scales.  Figure  14  is 
for  comparison  with  measured  temperatures  in  HEMSS  No,  3.  Figure  15  is 
for  comparisons  with  runs  made  on  the  3.66- and  6.10-m  (12-  and  20-ft) 
cylinders.  Results  in  Figure  14  compare  favorably  with  measured  tempera¬ 
tures  (Figure  11)  with  the  following  exceptions.  Comparisons  of  tempera¬ 
tures  at  gage  T1  should  not  be  made  since  the  temperatures  measured  by 
the  gage  are  assumed  to  be  air  temperature  while  the  calculated  values 
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at  the  same  position  represent  grout  surface  temperature.  Gage  T2  should 
be  used  to  assess  the  measured  and  calculated  values  near  the  inner 
surface.  Measured  and  calculated  temperatures  for  gage  T6  are  cc.nparable. 
The  measured  values  from  gage  T5  (Figure  11)  appear  to  be  lower  than 
should  be  expected  within  the  first  several  days  after  placement.  The 
temperature  should  be  much  closer  to  that  measured  by  gage  T4  since  the 
locations  are  symmetric  about  the  midpoint  between  surfaces  and  cooling 
conditions  are  comparable  at  each  surface  at  that  time.  The  calibrations 
at  higher  temperatures  for  T5  should  be  checked  if  possible. 

22,  Runs  21  and  22  were  made  with  the  3.66-  and  6.10-m  (12-  and 
20-ft)  diameter  cylinders,  respectively.  Peak  temperatures  reached  109  C 
(229  F)  and  116  C  (240  F) ,  respectively.  The  maximum  temperature  achieved 
at  a  placement  temperature  of  25.5  C  (78  F)  should  be  116  C  (240  F)  pro¬ 
vided  that  the  adiabatic  temperature  rise  curve  used  is  correct  after 
7  days  age.  Figures  16  and  17  show  temperature  histories  in  runs  21  and 
22.  In  order  to  show  the  effects  of  Increased  diameter  more  effectively 
temperatures  are  plotted  for  nodes  in  runs  21  and  22  that  are  located  at 
proportionally  the  scaled  distances  apart  as  all  the  gage  locations  in 
HEMSS  No.  3  and  run  20.  These  figures  show  the  increasingly  longer  time 
required  to  cool.  Also  evident  is  the  proportionally  larger  volume  cf 
grout  subjected  to  long-term  high  temperatures.  In  fact  peak  tempera¬ 
tures  in  the  6.10-m  (20-ft)  cylinder  have  only  lowered  to  66  C  (151  F) 
at  28  days. 


Thermal  Stress-Strain  Simulation 


23.  Thermal  stress-strain  simulations  of  runs  20-22  were  conducted 
using  the  University  of  California  program  described  earlier.  This  pro¬ 
gram  uses  time-dependent  mechanical  properties  and  thermal  loads  input 
from  the  temperature  program  to  calculate  thermal  stress  and  strain  at 
times  subsequent  to  placement. 

24,  All  mechanical  properties  were  assumed  for  2C4  grout  in  the 
simulation.  Modulus  of  elasticity  provided  by  Merritt  CASES,  Inc.,  for 
age  14  days  was  18,205  MPa  (2.64  x  10^  psi) .  A  linear  function  for 
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modulus  of  elasticity  versus  the  log  of  time  passing  through  the  14-day 
value  was  assumed.  Figure  18  shows  the  assumed  relationship.  Linear 
coefficient  of  thermal  expansion  was  estimated  to  be  3.33  x  10~^/°C 
(6.0  X  10  /°F) .  This  value  may  be  low  for  grout,  but  no  corroborating 
data  could  be  located,  especially  for  early  ages,  for  a  better  estimate. 
Poisson's  ratio  used  was  0.28. 

25.  It  was  assumed  that  creep  is  high  in  grout  especially  at  the 
temperatures  encountered.  Since  no  data  were  located,  estimated  values 
double  those  of  concrete  were  used.  Creep  values  are  input  to  the 
program  as  relaxation  curves  used  in  McHenry's  equation.^’  ^ 

N  /  -m. (t  -  T) 

£  (a,  t,  T)  =  a  2  A.(T)ll  -  e  ^ 

i=l  ^  \ 


where 

e  =  creep  strain 
c 

t  =  time  after  placement 
T  =  age  at  loading 
a  =  stress 

N  =  2  gave  satisfactory  fit  with  the  assumed  data.  Figure  19  shows 
relaxation  curves  for  creep  coefficients  A1  and  A2.  Constants  chosen 
for  ml  and  m2  were  0.34  and  0.52,  respectively. 

26.  When  creep  is  considered  In  the  program,  stresses  at  each 
time  step  in  the  analysis  are  modified  for  stress  relaxation  allowing  no 
strain  for  the  interval  up  to  the  next  analysis  time.  The  creep  parameters 
are  stored  and  the  change  in  stress  is  stored  as  residual  stress  to  be 
included  in  the  next  time  step  analysis.  When  these  stored  values  are 
applied  during  the  next  time  step  analysis,  strains  due  to  creep  are 
eliminated.  Thus,  strains  output  at  any  time  step  are  elastic  strains. 
Nodal  displacements  are  accumulated  and  reflect  the  total  distortion 
within  the  model.  This  requires  a  word  of  caution:  displacement  and 

not  strain  values  should  be  used  when  comparing  measured  versus  calculated 
distortion.  Elastic  strain  values  from  this  program  are  compared  against 
tensile  strain  capacity  to  assess  potential  for  cracking. 

27.  The  final  input  value  is  the  initial  stress-free  temperature 


16 


distribution  in  the  model.  This  Is  the  temperature  calculated  individually 
for  each  element  and  is  the  reference  for  computing  stress.  The  time  at 
which  stress-free  temperature  is  computed  is  che  age  where  elastic  proper¬ 
ties  begin  to  develop.  This  was  assumed  to  be  at  6  hr  for  2C4  grout. 

These  data  are  calculated  by  the  temperature  program  for  subsequent  use 
in  the  stress  program. 

28.  Four  thermal  stress-strain  analysis  computer  runs  were  made 
for  this  study.  Temperature  runs  20-22  for  the  2.13-,  3.66-,  and  6.10-m 
(7-,  12-,  and  20-ft)  diameter  cylinders  provided  thermal  loads.  Run  S20A 
was  an  analysis  of  the  2.13-m  (7-ft)  cylinder  without  creep,  and  runs 
S20-S22  were  analyses  of  all  three  cylinders  with  creep.  The  initial 

7  days  of  simulation  in  each  temperature  run  used  a  diurnal  temperature 
cycle,  and  thus  surface  temperatures  varied  with  the  changing  ambient 
temperatures  during  each  day.  Thermal  stress  analyses  are  more  than  an 
order  of  magnitude  more  costly  to  run  compared  to  the  temperature  calcu¬ 
lations.  Therefore,  it  was  necessary  to  use  longer  time  steps  between 
successive  calculations  for  each  analysis.  Except  for  increments  of 
0.25  days  and  0.5  days  used  for  days  1  and  2,  respectively,  all  remaining 
time  Intervals  were  1  or  2  days.  Thus,  the  effect  of  diurnal  cycle  upon 
thermal  stress  was  not  computed.  The  effect  calculated  was  approximately 
that  associated  with  mean  daily  surface  temperatures. 

2,13-m  (7-ft)  cylinder 

29.  The  results  of  runs  £20  and  S20A  are  the  thermal  stresses  and 
strains  presumed  to  exist  in  HEMSS  No.  3.  These  runs  as  well  as  S21-S22 
indicated  that  major  and  minor  principal  stresses  and  strains  were  so 
oriented  as  to  be  tangential  and  radial  stress  and  strain.  Stresses  and 
strains  produced  from  each  set  of  radial  rows  of  elements  were  identical. 
Since  comparison  with  estimated  tensile  strain  capacity  was  used  to  assess 
potential  cracking  and  since  internal  measurement  of  displf'cement  rather 
than  stress  is  made,  the  discussion  of  these  runs  is  based  upon  calculated 
strains. 

30.  Figure  19  shows  radial  temperature  gradients  and  the  resulting 
tangential  strain  gradients  for  analyses  with  and  without  creep  at  succes¬ 
sive  ages.  In  these  runs,  as  well  as  in  those  for  the  larger  cylinders, 


initial  peak  strains  occur  between  0.75  and  1.25  days  age,  The  strain 
distribution  is  characterized  by  highest  tensile  strain  at  the  inner  sur¬ 
face,  somewhat  lower  tensile  strain  at  the  outer  surface,  and  peak  com¬ 
pressive  strain  in  between.  Increasing  restraining  volume  with  increasing 
radius  contributes  to  the  higher  tensile  strains  at  the  inner  surface. 

In  both  cases  peak  tensile  strain  occurred  at  0.5  days  age  or  the  first 
change  of  temperature  relative  to  the  stress-free  temperature  at  0.25 
days.  It  is  also  evident  in  Figure  20  that  highest  surface  tensile 
strains  are  a  function  of  Intensity  or  slope  of  the  thermal  gradient  and 
not  the  total  temperature  differential.  Compressive  strains  peaked 
within  a  day  later  at  a  time  where  maximum  curvature  of  the  temperature 
gradient  at  peak  temperature  coincided  with  maximum  temperature  differ¬ 
ential.  Peak  tensile  strains  were  195  millionths  and  peak  compressive 
strains  were  104  millionths  with  and  without  creep.  As  time  passes 
strains  begin  to  decrease  due  to  stress  relaxation  and  cooling.  The 
largest  differences  with  and  without  creep  occur  near  temperature 
equilibrium.  For  example,  compressive  and  tensile  strain  at  9  days  with¬ 
out  creep  were  123  and  86  millionths,  respectively.  With  creep  these 
strains  were  65  and  46  millionths,  respectively. 

31.  Some  questions  remain  concerning  the  accuracy  of  the  creep 
correction  modeled  in  this  program.  It  appears  that  larger  stress  relaxa¬ 
tion  should  occur  at  very  early  age  especially  when  stresses  and  strains 
are  nearly  equivalent  to  yield  values.  It  is  recognized  that  additional 
research  and  computer  program  development  are  needed  in  this  area.  From 
the  beginning  of  the  study  it  was  agreed  by  all  parties  concerned  that  due 
to  lack  of  actual  mechanical  properties  and  due  to  known  computer  program 
deficiencies  these  results  should  not  be  expected  to  indicate  too  much 
more  than  a  comparative  pattern  of  behavior. 

3.66-  and  6.10-m 

(12-  and  20-ft)  cylinders 

32.  Runs  S21  and  S22  were  conducted  on  3.66-  and  6.10-m  (12-  and 
20-ft)  diameter  cylinders,  respectively.  Tangential  tensile  strains  at 
the  inner  surface  were  267  and  296  millionths,  respectively,  occurring 
at  0.75  days.  Figures  21  and  22  show  temperature  and  tangential  and 
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radial  strain  gradients  for  the  3.66-  and  6.10-m  (12-  and  20-ft)  cylin¬ 
ders,  respectively.  As  noted  earlier,  temperatures  in  both  cylinders  did 
not  reach  equilibrium  at  the  end  of  the  runs.  Therefore,  tangential 
stress  and  strain  are  continuing  to  Increase  in  compression  at  the  inner 
and  outer  surfaces  and  in  tension  in  between  as  cooling  continues.  As 
each  increase  in  diameter  produced  higher  strains  due  to  the  initial 
temperature  rise,  strains  produced  during  cooling  are  also  higher.  Radial 
strains  are  lower  in  magnitude,  but  respond  similarly  to  tangential 
strains  during  cooling. 

33.  Figure  23  shows  estimated  tensile  strain  capacity  used  in 
thermal  stress-strain  runs.  These  values  are  twice  those  measured  for 
mass  concrete.  These  values  may  be  substantially  lower  for  grout  but 
no  data  exist  to  establish  a  better  estimate.  Tensile  strains  plotted 
in  Figures  19-22  can  be  compared  to  values  for  tensile  strain  capacity 
in  Figure  23  at  equal  ages.  In  all  cases  surface  tensile  strains 
occurring  early  are  five  to  six  times  higher  than  estimated  strain 
capacity.  Tensile  strains  that  occur  within  the  mass  during  cooling 
are  all  lower  than  strain  capacity  at  the  corresponding  ages.  It  should 
be  stressed  that  the  analysis  is  elastic  which  may  contribute  to  sub¬ 
stantial  deviation  from  actual  conditions.  These  simulations  do  show 
comparative  effects  in  the  three  cylinder  sizes. 

Conclusions  and  Recommendations 

34.  As  a  result  of  this  study,  many  observations,  conclusions, 
and  recommendat ions  based  upon  both  analysis  of  results  and  review  of 
pertinent  literature  are  possible.  Two  basic  considerations  result. 

The  initial  concern  is  potential  excessive  thermal  stress  and  strain 
caused  by  the  restraining  mass  and  thermal  gradient.  The  second  is  the 
potentially  detrimental  or  unknown  effects  on  the  grout  due  to  the 
excessive  temperatures  during  early  curing. 

35.  This  study  showed  that  the  finite-element  temperature  pro¬ 
gram  in  use  at  WES  could  provide  sufficiently  accurate  results  of  tem¬ 
perature  rise  and  decay  in  the  HEMSS  cylinders  based  upon  assumed 
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thermal  properties  and  boundary  conditions  provided  that  the  time  rate 
of  heat  generation  or  measured  temperature  rise  are  available.  Evalua¬ 
tion  of  the  effectiveness  of  control  measures  can  be  made. 


36.  This  study  showed  that  peak  temperature  in  cylinders  at 
least  2.16  m  (7  ft)  in  diameter  can  only  be  effectively  controlled  by 
limiting  placement  temperature.  The  rate  of  heat  generation  is  so  high 
that  boundary  conditions  can  only  affect  the  rate  of  later  cooling  and 
cannot  control  peak  temperatures. 

37.  Calculations  of  thermal  stress  and  strain  based  upon  assumed 
mechanical  properties  for  2C4  grout  and  the  known  thermal  gradients 
existing  in  the  2 . 16-m  (7-ft)  HEMSS  cylinders  indicate  that  high 
tangential  tensile  strains  occur  at  the  inner  and  outer  surfaces  within 
the  first  day  after  placement.  Unless  the  actual  cylinders  show  other¬ 
wise  it  must  be  assumed  that  thermal  cracking  is  probable  at  these  loca¬ 
tions.  Peak  initial  stress  and  strain  occur  within  0.5  days  of  peak 
temperature  after  which  both  are  reduced  as  cooling  occurs.  Before 
temperature  equilibrium,  thermal  stress  and  strain  reach  zero  after 
which  stress  and  strain  increase  opposite  in  sign  to  that  occurring 
earlier  until  thermal  equilibrium  is  reached.  Areas  initially  in  tension 
subsequently  experience  compression  and  vice  versa.  Although  it  cannot 
be  confirmed  here,  conditions  exist  to  cause  early-age  surface  cracking 
to  propagate  into  or  through  the  radius  of  the  cylinder  as  cooling 
progresses  and  tangential  tension  builds  in  the  central  radius  area. 

38.  Thermal  simulations  of  3.66-  and  6.10-m  (12-  and  20-ft) 
cylinders  subjected  to  the  same  conditions  as  HEMSS  No.  3  showed  an 
increase  in  peak  temperature  reaching  a  probable  maximum  value  of  116  C 
(240  F)  in  a  6.10-m  (20-ft)  cylinder.  Thermal  stress  and  strain  continue 
to  increase  with  increasing  diameter. 

39.  Based  upon  this  study,  thermal  simulations  under  the  condi¬ 
tions  that  existed  in  HEMSS  No.  1-3  would  not  have  produced  satisfactory 
results  without  actual  data  from  measurements  of  temperature  rise  or 
adiabatic  heat  of  hydration  tests.  The  rate  of  initial  heat  of  hydration 
is  so  much  higher  than  heat  of  hydration  values  obtained  from  conventional 
tests  that  calculation  of  peak  temperature  would  have  been  at  least  22  C 
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(40  F)  below  measured  values.  In  addition,  expected  large  variations  in 
other  thermal  properties  assumed  constant  in  this  study  are  probably 
based  upon  concrete  data. 

40.  When  comparing  the  effect  of  creep  on  thermal  strain  calcula¬ 
tions,  the  results  showed  that  creep  does  not  affect  early-age  high  stress 
and  strain,  but  does  lower  the  final  stress  and  strain  at  thermal  equilib¬ 
rium  by  one-half.  This  is  true  only  if  the  assumed  values  of  mechanical 
properties  are  fairly  accurate  and  if  the  method  of  modeling  creep  is 
correct . 

41.  Some  investigators  believe,  based  upon  data  obtained  from 
measurements  of  strain  in  mass  concrete  structures,  that  stresses  at 
early  age  are  much  lower  than  should  exist  based  upon  elastic  properties 
and  geometry.  These  measurements  seem  to  show  that  as  soon  as  cooling 
begins  areas  supposedly  in  compression  begin  to  shrink  as  though  no 
compression  existed.  Concrete  that  is  restrained  from  expansion  does 
not  increase  in  compression  as  its  elastic  properties  would  indicate  at 
early  age.  The  result  is  somewhat  like  100  percent  creep  or  stress 
relaxation  during  early  stages.  Were  this  the  case  in  the  HEMSS  cylinders, 
tensile  and  compressive  strains  would  be  near  zero  as  soon  as  cooling 
began.  The  surface  would  immediately  undergo  compression,  the  interior 
mass  undergo  tension,  and  the  magnitude  of  these  values  at  thermal 
equilibrium  would  be  substantially  higher  than  calculated.  It  is  evident 
that  very  little  is  known  about  some  mechanisms  occurring  in  mass  place¬ 
ments  containing  hydrating  Portland  cement  at  early  age,  especially  at 
high  temperatures.  Therefore,  the  results  reported  here  may  be  subject 

to  substantial  error  under  the  circumstances. 

42.  The  second  major  concern  resulting  from  the  study  is  that 
2C4  grout  subjected  to  extremely  high  temperatures  at  early  age  may 
devel..p  substantially  different  properties  than  when  cured  at  ambient 
temperatures.  It  seems  necessary,  especially  if  larger  diameter  cylinders 
are  seriously  considered,  that  at  least  a  minimal  examination  of  mechani¬ 
cal  properties  pertinent  to  the  explosive  tests  be  obtained  for  2C4  cured 
through  the  range  of  temperatures  expected.  It  is  known  that  ultimate 
strength  is  reduced  as  placement  and  curing  temperatures  increase.  It 
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is  not  known  how  other  mechanical  properties  are  affected  at  high  initial 
curing  temperatures. 

A3.  If  minimal  testing  proves  that  pertinent  properties  change 
substantially  from  those  desired,  a  complete  test  program  to  evaluate 
2C4  as  a  function  of  temperature  is  warranted.  It  may  be  preferred  in 
that  case  to  choose  another  grout  with  lower  quantities  of  cement  and 
document  it  at  the  temperatures  expected.  Even  if  2CA  is  fully  documented 
the  thermal  stress  effects  from  high  temperature  gradient  are  still 
present . 

AA .  Questions  may  arise  concerning  the  advantages  of  adding 
pozzolans  such  as  fly  ash,  etc.,  to  replace  some  of  the  cement  in  2CA. 
Pozzolan  replacement  will  reduce  the  rate  and  amount  of  heat  of  hydration. 
Temperature  rise  may  be  reduced  by  20  percent.  Data  will  have  to  be 
obtained  for  the  revised  mixture.  Experience  at  WES  has  shown  that  pozzo¬ 
lans  become  less  effective  in  reducing  temperature  rise  when  the  cement 
content  is  high.  It  seems  as  though  another  mixture  with  lower  cement 
content  would  be  the  best  alternative  if  2CA  as  presently  designed  is 
unsatisfactory. 

A5.  No  substantial  decrease  in  peak  temperatures  can  be  obtained 
when  placing  2CA  grout  in  mass  unless  crushed  ice  replaces  the  mixing 
water  and  materials  are  precooled  to  result  in  a  placement  temperature  at 
or  slightly  higher  than  A. A  C  (AO  F).  Placement  at  low  temperatures  will 
retard  the  early.  Intense  rate  of  hydration.  Peak  temperatures  may  be 
lowered  more  than  the  difference  in  placement  temperatures  because  a 
lower  rate  of  hydration  will  cause  additional  time  to  occur  between 
placement  and  peak  temperature  allowing  additional  cooling  to  the  air. 

This  will  probably  only  be  possible  with  the  2.13-m  (7-ft)  cylinders. 

If  a  6.10-m  (20-ft)  cylinder  is  cast  at  the  placement  temperatures  used 
in  HEMSS  No.  3,  25.5  C  (78  F) ,  peak  temperature  will  be  approximately 
13.9  C  (25  F)  higher  than  in  HEMSS  No.  3.  Therefore,  reducing  the 
placement  temperature  to  A. A  C  (AO  F)  will  lower  peak  temperatures  only 
slightly  below  those  seen  in  HEMSS  No.  3.  Due  to  the  larger  diameter 
thermal  stress  and  strain  will  be  higher  than  in  HEMSS  No.  3  because 
volumetric  restraint  is  increased. 
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46.  Many  unanswered  questions  exist  on  the  use  of  2C4  grout 
placed  in  mass.  As  discussed,  the  questions  deal  with  the  material's 
temperature  dependence  and  quality  as  well  as  with  potential  excessive 
thermal  stress  and  strain.  To  answer  both  questions,  additional  thermal 
and  mechanical  properties  data  are  necessary  tor  2C4  grout  cured  at 
elevated  temperatures. 
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Figure  1.  Finite-element  model  No.  1  used  for  temperature  runs  1-19. 


FIGURE  3  -  ASSUMED  ADIABATIC  TEMPERATURE  RISE  FOR  2C4  GROUT 


FIGURE  7  -  DIURNAL  AIR  TENPERATURE  CVCLE  USED  IN  RUNS  10-12 
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Figure  13,  Radial  finite-element  node  locations  for  runs  20-22.  Location 
at  which  thermal  stresses  and  strains  are  calculated  for  elements  are 
radially  midway  between  adjacent  nodes. 
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Figure  21.  Computed  strain  and  temperature  gradients  across  the 
radius  of  3.66-m  (12-ft)  diameter  cylinder  at  successive 
times  after  placement  (Sheet  1  of  3) 
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PIGURE  23  -  f^SSUHED  TENSILE  STRAIN  CAPACITV  FOR  204  GROUT 
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